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GUIDING  CENTER  HAMILTONIAN  THEORY 
Of  ErEE-^Cectron  laseks 


Cel  so  Grebogi*  and  Han  S.  Uhtn 
Naval  Surface  Weapons  Center 
White  Oak,  Silver  Spring,  Maryland  20910 


ABSTRACT 


The  relativistic  guiding  center  ponderomotive  Hamiltonian  for  f ree¬ 
el  ectron  lasers  is  derived.  The  derivation  takes  into  account  arbitrary 
signal  wave  polarization  and  wiggler  field  geometry  including  guide  field 
nonuniformities.  In  particular,  it  is  allowed  (i)  for  a  tapered  axial  guide 
quasi  static  magnetic  field  along  its  direction  of  propagation,  (ii)  for  a 
realizable  wiggler  to  be  tapered  both  in  amplitude  and  in  period,  and  (iii) 
for  the  signal  electromagnetic  wave  to  be  a  growing  modulated  wave  of 
arbitrary  geometry  propagating  in  the  direction  of  the  static  magnetic 
field.  The  equations  of  motion  are  then  derived  including  the  guiding  center 
perpendicular  drifts  and  beam  quasi  static  self-fields. 


♦Permanent  address:  Laboratory  for  Plasma  and  Fusion  Energy  Studies 

University  of  Maryland 
College  Park,  Maryland  20742 
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I.  INTRODUCTION 


In  a  previous  paper^  (henceforth  referred  as  "I"),  a  general  expression 
for  the  relativistic  ponderomotive  Hamiltonian  of  two  interacting 
electromagnetic  waves  was  derived.  It  was  found  that  relativistic  effects 
introduce  new  terms  in  the  expression  for  the  ponderomotive  Hamiltonian 
similar  to  those  found  earlier^  in  the  case  of  a  single  electromagnetic 
wave.  The  derivation  in  “I"  was  as  general  as  possible.  It  was  allowed  for 
arbitrary  kp,  u/fl,  k^/kj, ,  v/c  <  1,  polarization,  and  for  slowly  growing 

and  spatially  modulated  waves.  It  was  necessary  to  use  Hamiltonian  theory  in 
order  to  account  for  all  that  and  to  introduce  nonuniformities  in  the 
background  fields  in  a  systematic  way.  The  calculation  was  an  application  of 
relativistic  guiding  center  theory^  modified  by  the  averaging  over  the  two 
interacting  electromagnetic  waves. 

The  theory  developed  in  "I"  is  suitable  for  free-elect ron  lasers^  with  a 

guiding  magnetic  field  in  which  the  gyrpperiod  is  comparable  to  the  period  of 

the  signal  wave  and  to  the  time  it  takes  for  the  electron  beam  to  travel  one 

wavelength  of  the  wiggler  [fl  ~  ui  ~  k^  operating  regime 

ii 

which  gives  maximum  efficiency.'*  The  guiding  center  theory  of  the  nonlinear 
response  of  a  relativistic  particle  beam  in  a  guide  magnetic  field  to  the 
signal  wave  and  the  wiggler  field  has  not  been  adequately  treated  in  its 
generality.  In  this  work,  we  reduce  the  general  relativistic  ponderomotive 
Hamiltonian  expression  obtained  in  "I"  to  free-electron  lasers  and  derive  the 
equations  of  motion. 
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We  should  mention  that  the  Hamiltonian  formulation  of  free-electron 
lasers  has  the  advantage  of  expressing  the  vector  evolution  equations  in  terms 
of  a  single  scalar  function  on  phase  space.  Furthermore,  since  the  guiding 
center  equations  have  been  averaged  over  gyrations,  they  are  much  more 
suitable  for  use  in  numerical  integrations  and  simulations  of  free-electron 
lasers.  Namely,  instead  of  integrating  along  the  gyrating  trajectories  for 
each  of  the  electrons  in  the  beam,  it  is  possible  now,  for  the  averaged 
equations,  to  take  a  two  orders  of  magnitude  larger  integration  steps  with 
equivalent  savings  in  computer  time.  In  addition,  because  the  guiding  center 
transformation  has  been  done  in  the  Hamiltonian,  the  derived  equations 
conserve  phase  volume  exactly  and  thus  enhancing  the  equations  predictability 
power.  Finally,  although  we  are  dealing  with  single  particle  motion,  the 
derivation  of  the  averaged  Hamiltonian  is  the  major  part  of  the  work  for  a 
derivation  of  the  nonlinear  relativistic  gyrokinetic  equation  for  two 
interacting  waves.  Once  we  have  the  ponderomoti ve  Hamiltonian  for  the  free- 
electron  laser  in  a  guide  magnetic  field,  the  Vlasov  equation  for  the 
distribution  of  guiding  centers  F  in  the  Poisson  bracket  form 

K|.0, 

is  the  gyrokinetic  equation.®  This  will  be  reported  in  the  future. 

In  Sec.  II,  we  present  a  summary  of  the  essential  results  derived  in 
'T'.  Because  the  formulas  presented  here  are  so  general  we  can  reduce  them  to 
arbitrary  signal  wave  polarization  and  any  wiggler  field  geometry  including 


guide  field  nonuniformities.  Hence,  in  Sec.  Ill,  as  an  illustrative 
application  of  our  theory,  we  obtain  the  guiding  center  relativistic 
ponderomotive  Hamiltonian  for  free-elect ron  lasers  by  considering: 

(i)  the  axial  guide  quasi  static  magnetic  field  to  be  tapered®  along  its 
direction  of  propagation; 

(ii)  the  realizable^  wiggler  field  to  be  tapered®  both  in  period  and 
amplitude; 

(iii)  the  signal  electromagnetic  wave  to  be  a  growing  modulated  wave  of 
arbitrary  geometry  propagating  in  the  direction  of  the  static  guide 
field.® 

In  Sec.  IV.,  we  derive  the  guiding  center  free-electron  laser  equations 
of  motion  for  the  Hamiltonian  derived  in  Sec.  III.  including  the  guiding 
center  perpendicular  drifts  and  beam  quasi  static  self-fields.  A  conclusion  is 
presented  in  Sec.  V. 
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I I .  RELATIVISTIC  PONDEROMOTIVE  HAMILTONIAN 
Cf  two  interacting  waves 


As  derived  in  "I",  the  two  electromagnetic  waves  are  represented  by  their 
scalar  and  vector  potentials  (|)  and  A  in  an  arbitrary  gauge.  They  are  assumed 
to  have  the  form 


(t.(x,t) 


2  i'l'.(x,t) 

I  ?.(x,t)e  "  +  C.C.. 

i=l  ^ 


(la) 


and 


2  i'l>.(x,t) 

A(x,t)  =  A.(x.t)e  "  +  c.c. 


(lb) 


where  the  over-tilde  represents  the  slowly  varying  amplitude  of  the  wave 
packet.  The  local  wavenumbers  kj  and  ^2  and  frequencies  and  of  the 
waves  are  given  by 


k.(x,t)  =  Vi|».(x,t),  i  =  1,2. 


(2a) 


and 


(x,t) 

I.  [x,t)  =  ,  i  -  1,2 


The  relativistic  particle  Hamiltonian  is  given  by 


(2b) 
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(3) 


H(x,g)  =  {m^c^  +  c^[g  -  I  Aj^(x,t)  -  I  A[x,t)]}^2 
+  e<|)j^(x,t)  +  e(t>(x,t), 

where  we  allow  the  unperturbed  fields,  as  given  by  and  Aq,  to  be  slow 
functions  of  position  and  time.  Physically,  this  means  that  the  scale  length 
of  the  spatial  inhomogeneities  in  (j)^  and  Ag  (or  Eg  and  Bg)  are  smaller  than 
the  gyroradius,  which  is  the  standard  approximation  in  guiding  center  theory 
and  true  in  free-electron  lasers.  In  addition,  it  means  that  we  allow 
(j>g  and  Ag  to  change  appreciably  on  a  drift  time  scale. 

Because  we  are  interested  in  the  pond.eromotive  Hamiltonian,  which  is  a 
second  order  effect,  we  expand  (3)  up  to  second  order  in  the  waves 
amplitudes.  The  result  is 


(4a) 

(4b) 

(4c) 
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As  discussed  in  Ref.  2,  we  introduced  a  velocity-like  variable  Uq  by 

"lio  °  e  ■  f  5o "  ”ii  *  1 5* 

where  u  =  yv  is  the  world  velocity  and 


YqIiic^  is  the  energy  of  the  injected  relativistic  beam.  The  introduction  of 
Uq  allows  us  to  exhibit  the  perturbation  in  the  Hamiltonian  explicitly  instead 
of  in  the  Poisson  brackets. 

Next  we  transform  the  Hamiltonian  (4)  to  the  relativistic  guiding  center 
variables,  by  the  transformation  (x,Uq)  (X,U||,u,0).  The  first  three 
variables  X  represent  the  location,  in  physical  space,  of  the  particle's 
guiding  center.  The  variable  U|,  is  called  the  parallel  world  velocity  of  the 
guiding  center;  it  is,  except  for  the  relativistic  factor,  the  parallel 
guiding  center  drift  to  lowest  order  in  the  waves  fields.  The  next  variable 
u  is  the  magnetic  moment  and  is  related  to  the  perpendicular  guiding  center 
drift  by  u  =  muJ/2Bg,  where  Uj^  is  essentially  the  gyrophase  average  of  Ug^. 
Finally,  0  is  the  gyroaveraged  gyrophase  whose  time-derivative  to  lowest  order 
is  the  local  relativistic  gyrofrequency.  The  new  guiding  center  variables  are 
defined  by^ 

mcUg 

X  =  X  -  e  — a  +  O(e^)  ,  {5a) 

-  -  eog 

U,.  =  +  0(e),  (5b) 

Uj 
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+  0(e). 


(5c) 


and 


mur 


2B, 


0  =  0  +  0(e), 


(5d) 


where  the  fields  are  evaluated  at  (x.t)  and  e  is  the  dimensionless  parameter  t 

which  indicates  the  order  of  various  terms  in  the  guiding  center  expansion. 

Since  in  this  section  we  are  presenting  only  results,  we  will  set  e  =  1  from 
now  on  to  recover  the  physical  formulas.  The  instantaneous  gyrophase  0  is 
defined  implicitly  by 


a  =  cos  exj  -  sinex^,  (6a) 


and 


A  A 

c  =  -  sinSx^  -  cos9x2.  (6b) 

A  A  A 

The  perpendicular  unit  vectors  a  and  c  rotate  with  the  particle,  a  is  in  the 

A  A  A 

direction  of  the  gyroadius  as  given  by  (5a)  and  c  is  given  by  Uq  =  Ug  b  +  Uq  c, 

A  AAA  IIX 

where  b  =  Bg/Bg.  The  triad  (a,b,c)  obeys  a  =  b  x  c.  The  perpendicular  unit 

AAA 

vectors  do  not  rotate  with  the  particle  and  satisfy  x^  x  Tj,  =  b.  They  have  a 
slow  dependence  on  (x,t)  since  Bg  itself  does. 

% 

For  the  Hamiltonian  structure  to  be  complete  we  need  the  Poisson  brackets  4 

among  the  relativistic  guiding  center  variables.  The  guiding  center  variables 
are  noncanonical  and  the  nonvanishing  Poisson  brackets  are  given  by^ 
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(7a) 


+  ^  fills  _  ii)  .  (9) 

me  '^30  3w  3u  30^  *  ' 

This  formula  will  be  useful  in  deriving  the  equations  of  motion  from  the 
guiding  center  Hamiltonian. 

In  terms  of  the  guiding  center  variables  (5),  the  unperturbed  Hamiltonian 
(4a)  becomes 
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(10) 


Hq  =  64)^  +  me  r , 


where 


r  ^  (1 


2yBrt(X,t)  U 


The  first  order  Hamiltonian  (4b)  is  given  by  the  following  Fourier  expansion^ 


2  +«  i4',.[X,t)+iil[0  +  a.(X,t)  + 

H,  =  I  I  {H.  e  ’  '  ■  ^  +  c.c.}, 

^  i=l  Jl=-“  l 


where  H.  is  the  Ath  9 -Fourier  coefficient  of  the  first  order  perturbation  due 

%  — 
to  the  ith  wave  and  a.(X,t)  is  the  angle  between  and  .  We  introduce  the 

triads  (k.  ,  b  x  k.  ,  b)  for  i  =  1,  2,  in  which  k.  *  cosa.  +  sin  a. 

U  'i  -L 

The  Fourier  coefficients  are  then  given  by 

j  ji 

^i  " I ?i*  ^  3ir  ^ 


where  j'[  a  J^(k.  p)  ,  p  =  \  a  =  and  7,^  =  U„/r.  As  for  the  second 

order  Hamiltonian,  we  only  ke8p  the  slow  terms  in  Eq.  (4c)  since  we  are 
interested  in  the  ponderomotive  Hamiltonian.  We  average  over  gyration  and 
oscillation  of  both  waves  but  keep  the  slow  beating  terms  of  the  form 


expi['|»2(X,t)  -  4'2(X,t)].  The  result  is^ 


lil  -42' 

r  me 


i('pT  -  'l'2) 


+  C.C.]  (14) 


+  -7  KpSq  +  [|A,  I^IAj  f  t  Jg{Akj^p)(Aj  •  Aj 


i('l'T  -  '1^2- 


i  “i 


■5  9  ~  ~  *  U'h-’l'o' 

ITT  i2,nr  Zj.  1  tai,  a  ^  '  2 


+  c.c.)  +  2vBq  [|A^  I^+Ia^  P+  jQ(Ak^p)(A^  Ag  e  '  ^  +  c.c.)]}). 
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where  Akj^  =  [k^  +>2  -  2k-|  k^  cos[aj  -  a^)]  and  we  made  use  of  the 

X  «!•  X  J» 

Bessel  identity 


At  this  point  we  subject  the  Hamiltonian  H  =  Hg  +  to  an  averaging 

transformation  to  eliminate  gyration  and  fast  oscillations  contained  in  but 
keeping  the  beating  terms  expi(i|^j  -  This  canonical  transformation  is  effected 
by  using  Lie  transforms  which  leads  to  the  following  0-independent  Hamiltonian 


K  =  Kg  +  K2. 

where  Kg  =  Hg  as  given  by  Eq.  (10),  or 
Kg  =  e(j>g  +  mc^r. 


(15) 


(16) 


and  the  ponderomotive  Hamiltonian  K2  contains  two  parts 


Kj  -  +  K^. 


The  first  part  K^  =  H^  which  is  given  by  Eq,  (14).  The  second  part  comes  from 
the  canonical  transformation  of  Hj^  and  can  be  expressed  as^ 


2  +  »  I  ^2  ®  '  '  +  c,c.] 

y  y  r®i  i_  +  —  J—]  ^  ^  ^ _ 

^mc  3u  m  3U.|  0).  -  k.  M  -  IQ. 

1  •*  I  Xk****®®  I!  1  <  Jl  ” 
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Where 
expii(a^ 
waves  is 


and  are  given  by  Eq.  (13).  It  is  possible  to  absorb 
i  * 

-  a'}  into  .  In  this  case  the  triad  of  unit  vectors  for  both 

A  A  A  ^  A 

(k.  ,  b  X  k,  ,  b)  in  which  we  drop  the  subscript  I. 


V 
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III.  RELATIVISTIC  PONDEROMOTIVE  HAMILTONIAN 
fOr  free-electrqn  lAse^ 


We  consider  a  relativistic  electron  beam  propagating  in  the  tapered  guide 
field  Bq  =  8Q(Z,t)b,  The  free-electron  laser  effect  results  from  the  bunching 
of  the  beam  due  to  the  beating  of  the  signal  wave,  represented  by  the 
potentials  <j)  and  A  ,  and  the  magnetostatic  wiggler  A  .  For  the  purpose  of 
illustrating  our  theory,  we  reduce  the  relativistic  guiding  center  Hamiltonian 
of  two  interacting  waves,  as  given  by  Eqs.  (13)  -  (18),  by  assuming  the 
following  signal  and  wiggler  fields  geometry: 

(1)  The  signal  electromagnetic  wave  is  a  growing  modulated  wave  of 
arbitrary  geometry  propagating  in  the  Z  direction, 

i/  [ks  (Z')+  "k*  (Z')ld2'  -  1»t 

('JX.t)  =  I  I,  (X,t)e  Oil  .  c.c.  (19) 

"  n  ” 

and 

Z 

i/  [kj  (z'l  +  nk^  (Z')]dZ'  -  1<ot 

♦JX.t)  ■  I  ?,(x,t)e  Oil  ^  ,21,, 

^  “  n  ^  ■ 

where  the  sum  is  over  the  wiggler  harmonics  and  the  subscripts  1  and  ii 

A 

indicates  perpendicular  and  parallel  to  the  local  b,  respectively. 


(2)  The  realizable  wiggler  field  is  tapered  both  in  period  and  amplitude 


+  c.c 


(21) 


i/\  (Z‘]dZ' 

As  indicated  by  Eq.  (17),  we  split  the  ponderomoti ve  Hamiltonian  in  two 
parts.  The  first  part  is  given  by  Eq.  (14)  and  its  reduction  for  free 
electron  lasers  is  straightforward  by  identifying,  say,  =  A^  and  ^  =  A^. 
Then  for  the  signal  and  wiggler  fields,  as  given  by  Eqs.  (19)  and  (21), 
jQ(Ak^p)  =  1  and  we  obtain 


+  2  I  |A  •  a*  I  cos  (/  [k'(Z')  -  k  (Z')]dZ'  -  ujt) ,  (22) 

n  0  ’^11 

where  r  is  given  by  Eq.  (11)  and  k*  =  k  +  nk  .  We  observe  that  the  fields 

«  S||  W|| 

amplitudes  Ja^  |  and  |a^  |  depend  on  (X,t). 

The  second  part  is  given  by  Eq.  (18)  where  H  and  H  is  given  by 
Eq.  (13).  The  reduction  of  for  free  electron  laser  fields  (19)  -  (21)  is 
more  involved  than  for  K^.  For  both  signal  and  wiggler,  only  the  A  =  0  and 
A  =  +  1  Fourier  components  are  different  from  zero,  as  given  by 


-  0  “ 

(23a) 

C 

H  1 

0, 

(23b) 

V 

o 

It 

Hs 

^  A  *  ±  1 

(23c) 
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and 


z  =  ±  1 


1  e  -  Tf  _  in  Tf 


(23d) 


where  we  absorbed  expiil(a  -  o  ]  in  H  ,  e  is  the  signed  electron  charge,  and 

s  w  w. 

A  A  AAA  ^ 

we  chose  x  =  k,  and  y  =  b  x  k,.  Therefore,  A  ,  A_  and  A  ,  A  are  in 

^  lx  ^ 

relation  to  the  x,y  frame.  In  view  of  expressions  (23),  K2  can  be  written  as 
the  Cartesian  components  of  the  vector  potentials 


kI  ■ 


4=0 


lA  =  1 


lA  =  -  1 


which  are  obtained  from  Eq,  (18).  For  A  =  0,  we  obtain 


In  taking  the  derivative  in  Eq.  (25)  we  should  keep  in  mind  that  V||  =  U||/r 


where  r  depends  on  U,  as  given  by  Eq.  (11).  Thus,  we  obtain 


1  =7  -  ^ 

=  0  n  rm(a)  -  k.lvj'^ 


To  evaluate  K, I  ,  we  need  the  absolute  values  of  Eqs.  (23c)  and  (23d)  . 

=  +  1 

and  their  cross  product  as  given  by 


lA  =  ±  1 


2  2 

4  ,2  l^sj  ^ 


(27a) 


2  2 

il^A  I' 

4  ^2  iwj  ^ 


(27b) 
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and 


H 


where  p  =  V 


K 


where 


R 


and 


R 


1)1  =  ±  1  '^U  =  ±  1 


/y  \ 

•  A  +  1  A, ,  J , 
-w^  -Sj. 


(27c) 


/ft.  Hence,  substituting  expressions  (27)  in  Eq.  (18)  we  obtain 


k! 


1  -  Y  r  r-i  i _ +  ^  —1  ^  _ 

=  ±  1  ~  n  “  -  k|;V|,+  ft 


e  3 


2  2 

i!AR* 

4  2  ^S 


e^V^ 

1  i  ,^+ 

*^w  ?  ^ 

_3_  +  l_  L  ‘ 

'•  m  3U||  “  nic  3ii" 


H  ^  w 

V.  +  ’ 


'w„  « 


(28) 


i  -=  |Asj'-^l?s^  *  -  k./Z')]  dZ-  -  .t) 

±  |A  A*  |sin[/  [k;(Z')  -  k  [Z')]dZ'  >  wt) ,  (29a) 

'-Si  -Wi'  Q  V»|| 

I  -  dz’  -  »t) 

t  -  t  I  S1"  (/[k|(Z’)  -  k  (Z')]  dZ’  -  a,t).  (29b) 

-Si  -Wi  Q  il  Wj 


V 
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To  evaluate  the  derivatives  in  Eq.  (28)  we  should  again  realize  that  Vj^  and  Q 
depend  on  r  which  in  turn  depends  on  u  and  U|| ,  Taking  that  into  account,  we 
obtain 


K«|  -  I  - ^  {[1  +  - 72  j  -'w 

=  ±  1  n  4rmc^  (-k  V.,  + 

W„  II 


• 

"  rl  R? 


-k;v,  *  !1  kjV  -  (-k;v,T 

*  U  <•  -  k;V,  i  n  .  k|v,  .  a)*'  • 


(30) 


Substituting  Eqs.  (26)  and  (30)  into  Eq.  (24),  we  obtain  the  expression  for 
the  second  part  of  the  ponderomoti ve  Hamiltonian  which  is  given  by 

i/b  _  Y  nic 

Kp  -  2,  .  ,,2 


2  2 


i^i[i  + 


k^ 

w,. 


n  4rmc 


[1 1 


k^ 


In-  To  ’**^“'^1' 

2^K  *  w  -  k‘V„ 


-  n 


0 


kj  c'  [-kJV||  -  £l)  . 

(-k;v„..)^ 


+  r?  "'^11  ^  +  kj^c-  -  9<;Vn  +  n)^  ^ 

“  ■  ^  (o)  -  k,;v„  +  fl)^  ^  * 


(31) 
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IV.  FREE-ELECTRON  LASER  EQUATIONS  OF  MOTION 


The  equations  of  motion  for  the  relativistic  guiding  center  Hamiltonian 
resulting  from  the  beating  of  the  wiggler  and  signal  fields  are  found  with  the 
help  of  Eq.  (9).  We  obtain  six  equations  of  motion  corresponding  to  each  one 
of  the  dynamical  variables. 


A,  Guiding  center  magnetic  moment. 

=  {vijK}  =  0,  or,  y  =  constant.  (32) 


Eq,  (32)  represents  the  adiabatic  invariance  of  the  magnetic  moment;  it  is  not 
an  exact  invariant  because  the  ponderomoti ve  Hamiltonian  still  depends  on  the 
gyrophase  to  order  cubic  in  the  fields.  As  it  was  shown  in  "I",  the  adiabatic 
magnetic  moment  can  be  expressed  in  terms  of  mUj^/2Bg,  which  contains  the  fast 
oscillations  of  the  signal  and  wiggler,  by 


y 


I 

=  ±  1 


^  H  exp(i$  ) 

r  \  \ 


+  I 

n 


I  Hg  exp(i$g  ) 


(33) 


where 


\  t  [0  +  a^[X.t)  +  Tr/2]  , 


and  similar  expression  for  $  . 

h 


i 


<1, 
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B.  Nonlinear  gyrofrequency. 


s-le.K) 


The  last  term  is  the  nonlinear  gyrofrequency  shift  due  to  the  ponderomoti ve 
potential.  The  w  derivative  of  Eq.  (22)  gives 


me 


r  2e^n 

i  r2„2.4 
n  r  m  c 


(1 


z 


(35) 


We  use  Eq,  (31)  to  evaluate  3K2/9y,  which  gives 


k'V„ 


1^2  ^2 


2  2  4 
n  2r^m^c^ 


k'^c^ 

vfn  y 


±)  R+  +  ( IL.  +  “ 

^  S  ^0) 


k'^c^  -  r^ 
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-]  r;} 


- 1 


^  {[3 


k2  c2 
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n  4r^m^c® 


k2  c2 

2I 
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8^ 

■  “+ 

+ 

3 

T 

8r 

7r^ 

2r^ 

Kj  1. 

0) 

0.2 

2r 


=  -  k,jV||  ?  and  0)^  =  (0  -  kjVn  +  n. 


The  sum  of  Eqs.  (35)  and  (36)  is  the  nonlinear  gyrofrequency  shift 
indicated  in  Eq.  (34). 

C.  Parallel  force  on  a  guiding  center. 

mU„  =  {mUj,  Kg  +  K2}  »  -  b  •  7(Kq  +  Kg).  (38) 

From  Eq.  (10)  we  obtain 

-  b  •  VKg  *  b  •  (eEg  -yVBg),  (39) 

where  Eg(X,t)  is  the  electric  field  due  to  the  beam  self-fields  and  Bg(Z,t)  is 
the  tapered  guide  field.  From  Eqs.  (22)  and  (31)  we  obtain  the  contribution 
of  the  ponderomotive  Hamiltonian  to  the  parallel  force,  which  is  given  by 


2  liB  Z 


n  rmc 
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and 


-  l<„  }  I 
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l<„  )dZ'  -  ‘-til. 
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_  + _ !! _ +  !!±  +  !L: 
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D.  Guiding  center  drifts 


(42) 


^  8K 

5=1^“  (Ti“''®0  -  '5n)  +  ^  (v,  ti  j/). 


II  m  311,. 


(43) 
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The  first  term  in  Eq.  (43)  are  the  perpendicular  drifts.  The  second  term  is 
the  guiding  center  parallel  velocity  which  differs  from  Vj  by  the 
ponderomoti ve  term  correction.  Explicitly,  (1/m)  is  evaluated  using 

expressions  (22)  and  (31)  which  yields 
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1  ^*^2 


m  3U 
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(44) 


i^^2  . 

i  w  ■  ‘  ^  t::  i  ^  .2  J  L  _2  „2 


2k,; 


II  n  2„4 


m'^r  (o)  -  k'V  ) 


c^  r‘^(a)  -  k;vj 


•]. 


3e^V?V 


k2  c2 
w,. 


n  4r^/c® 


?!  I'  * 


C*“ 

^w. 


2!  «;• 


2r 

+  0). 


ri  2r  kj^c^ 

*  C  *  -  V 


-  r 


0) 


-)  r;i. 


+  i 


"  ■(-  k^^v,  -  3)‘  -k^^v, 


(- 


2k2 

*1 


(■ 


k  V„ 
w,  II 


^  4) 


-  n  c 


2k2  c2 


k  V,, 
w„  II 


«V, 


r  " _  .  _±'\  D“ 

''-k..  V„  +  n  .2''  w’ 


w 


22 


(45) 
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2a)^(-k;v 
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The  sum  of  Eqs.  (44)  and  (45)  are  the  ponderomoti ve  correction  to  the  parallel 
velocity  drift. 


Equations  (32),  (34),  (38),  and  (43)  are  the  complete  guiding  center 
equations  for  free  electron  lasers.  Having  these  equations  at  hand,  we  can 
obtain  the  equation  for  the  energy  variation  of  the  particles  which  is  given 

byi»2 


dE  _  ^^*^0  *^2^  e  ^  .^^TOTAL 

W  ~  3t  "  c  -  3t  » 


(46) 


fn  ^  9 

where  A-j-Qjy^L  “  -O  -s  -w  ^  e  Qiven  by  Eq.  (43).  Similarly, 

we  can  straightforwardly  write  the  equation  for  the  phase  variation,  which  is 
gi ven  by 


+  -k„(X,)l  X,. 

11 

where  Xj  is  obtained  from  Eq.  (43). 


(47) 


23 


V.  CONCLUSIONS 


We  have  obtained  a  general  expression  for  the  ponderomotive  Hamiltonian 
for  free-electron  lasers  as  given  by  Eqs,  (22)  and  (31).  From  this  expression 
we  derived  the  general  guiding  center  equations  of  motion  governing  the  full 
free-electron  laser  phase  space  dynamics.  They  can  be  simplified  depending  on 
the  specific  experimental  system  under  consideration.  For  instance,  whole 
terms  will  drop  out  just  by  assuming  uj  «  uj.  Further  simplifications  can  be 
done  by  keeping  the  terms  with  the  right  resonant  denominator  for  a  given 
problem. 

We  should  also  mention  that  although  we  assume  X  to  be  Cartesian  in  our 
derivation,  we  can  express  X  in  terms  of  other  systems  of  coordinates,  say, 
cylindrical,  flux,  or  Hamada,  by  the  application  of  the  chain  rule. 

Our  final  expressions  are  given  in  terms  of  the  scalar  potential  <t>(X,t) 
and  the  perpendicular  vector  potentials  ^  (X,t)  and  It  Is  necessary 

to  choose  the  appropriate  form  of  these  realizable  potentials  for  a  given 
experimental  situation. 
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